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Abstract 

A computational study was conducted to better 
understand experimental results obtained from wind tunnel 
tests of a Mach 4 waverider model and a comparative refer- 
ence configuration. The experimental results showed that the 
performance of the reference configuration was slightly better 
than that of the waverider model. These results contradict 
waverider design theory, which suggests that a waverider 
optimized for maximum lift-to-drag should provide better 
performance than any other non-waverider configuration at a 
given design point, especially at hypersonic speeds. The 
computational results showed that the predicted surface pres- 
sure values and the integrated lift and drag coefficients from 
the pressure distributions were much lower for the reference 
model than for the flat-top model, due to the reference model 
bottom surface having a slight expansion. The lift-to-drag 
ratios for the flat-top model were higher due to a relatively 
low drag for the same amount of lift. These results indicate 
that the performance advantage of the reference model was 
due to the shape of the bottom surface and not due to the flat 
top surface. The results also showed that the reference model 
exhibited the same shock attachment characteristics as the 
waverider because the planform shapes were identical. CFD 
predictions show that the planform shape gives the waverider 
an advantage in performance over conventional hypersonic 
vehicles and that altering the bottom surface of a waverider 
does not cause significant performance degradation. 

Nomenclature 

a angle of attack, degrees 

Cl lift coefficient 

Clj drag coefficient 

X,Y,Z coordinate axis system, inches 
I.J.K computational grid axes 

M Mach number 

P static pressure, lbf/ft~ 

Pu ase base pressure, lbf/ft" 

\' inner law variable 

Subscripts 

°° freestream conditions 


* Aerospace Engineer, Supersonic/Hypersonics Aerodynam- 
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Introduction 

Waveriders are candidate shapes for various types of 
hypersonic aircraft designs. A waverider is a shape designed 
from a known supersonic or hypersonic flow field, such as 
flow past a right circular cone. 1 The waverider is designed 
such that the bow shock is attached along the outer leading 
edge at the design point. A typical waverider design from a 
conical flow field is shown in figure 1 . The lower surface is 
designed by using the known flow field to trace streamlines 
from the leading edge to the trailing edge. A waverider is 
uniquely defined by a leading edge definition and a specific 
set of freestream conditions. The upper surface may be 
designed as either a freestream surface or as a slight expan- 
sion surface to provide an additional contribution to lift. The 
shape is generally optimized for either maximum lift-to-drag 
ratio or minimum drag at the design point. The attached 
shock wave creates an efficient compression lifting surface 
with no flow spillage from the lower surface to the upper sur- 
face at the design point. Because of this characteristic, the 
predicted lift-to-drag ratios of waveriders are higher than 
those for conventional hypersonic vehicle concepts. The uni- 
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Figure 1 . Waverider Designed From Conical Flow Field 
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formity of the lower surface flow field and the absence of 
crossflow make the waverider an ideal candidate for scramjet 
propulsion system integration.^ 

Previous criticisms of waveriders have been 
resolved in recent studies, 13 4 leading to renewed interest in 
their use for hypersonic vehicle designs. One criticism is that 
early design techniques used only inviscid methods which 
produced shapes with large surface areas that resulted in large 
skin-friction drag. The lift-to-drag ratios obtained experimen- 
tally were poor compared to predictions. Current design 
codes include an estimate for skin-friction drag in the optimi- 
zation process. 1 This improvement allows the optimization 
routine to minimize wetted surface area and provide more 
accurate estimates for aerodynamic performance. The class 
of vehicles designed using this method is called viscous-opti- 
mized waveriders. A second concern associated with waver- 
iders is that their off-design performance may be poor even 
though on-design performance is excellent. However, recent 
experimental and computational studies have shown that 
waveriders which demonstrate acceptable off-design perfor- 
mance can be designed. 3,4 Additionally, optimization rou- 
tines can include various volumetric constraints which allow 
for the design of shapes with improved volumetric efficien- 
cies and packaging characteristics while accepting a mini- 
mum penalty in aerodynamic performance. Previous 
computational studies have validated design methods and 
predicted waverider flow-field properties accurately. 4 

An experimental program was conducted at NASA 
Langley Research Center to investigate the aerodynamic per- 
formance of a viscous-optimized waverider and a reference 
configuration. 3 A Mach 4 waverider and a flat-top reference 
configuration were tested in the Unitary Plan Wind Tunnel 
(UPWT). The design point of Mach 4 was chosen based on 
facility limitations and the desire to obtain data at, above and 
below the design Mach number. The reference configuration 
was intended to show the benefits of waverider flow-field 
properties by comparing a waverider to a non-waverider at 
the same conditions. The configurations were tested over an 
angle-of-attack range from -16° to 14° at Mach 4 and at off- 
design Mach numbers at zero angle of attack. This study 
showed that the maximum lift-to-drag ratio of the flat-top 
model was approximately 5 percent higher than that of the 
waverider model. This result was unexpected since theoreti- 
cally the waverider should provide better performance than 
any other non-waverider configuration at the design point. 

In order to better understand these experimental 
results, a computational study was conducted. The study con- 
sisted of obtaining computational fluid dynamics (CFD) solu- 
tions of both the waverider and the reference model flow 
fields. The results were used to interpret the wind tunnel data 
for the two models and to illustrate the benefits of waverider 


flow-field properties. This paper will present the details of 
the CFD study and limited comparisons with the experimen- 
tal data. 

Configuration Design 

The waverider model was designed using the Mary- 
land Axisymmetric Waverider Program (MAXWARP). 1,5 
The MAXWARP code includes an estimate for skin friction 
in the optimization process, using the reference temperature 
method. The code also uses a simplex optimization algorithm 
to optimize shapes for maximum lift-to-drag ratio or mini- 
mum drag at the design point. 5 In this study, the waverider 
model was optimized for maximum lift-to-drag ratio at Mach 
4.0. This is more appropriate as a hypersonic cruise perfor- 
mance parameter than minimum drag. The top surface was 
designed as an expansion surface using the axisymmetric 
method of characteristics. Some volumetric constraints were 
also incorporated into the optimization routine in order to 
increase volumetric efficiency and to generate a shape which 
had good structural characteristics for a wind tunnel model. 
These included a range for the base-height-to-length (fine- 
ness) ratio and a minimum volume. The cone semi-apex 
angle used to create the conical flow field was chosen to yield 
the highest maximum lift-to-drag ratio. An oblique view and 
a 3-view sketch of the waverider model is shown in figure 2. 

The flat-top model was designed by adjusting the 
top surface coordinate of the waverider at each cross section 
to create a flat surface and then adjusting the corresponding 
bottom surface coordinate at that cross section by the same 
increment. The result is a non-waverider configuration with 
the same cross-sectional area distribution and planform shape 
as the waverider model. An oblique view and 3-view sketch 
of the reference model is shown in figure 3. 

Computational Study 

The computational study consisted of obtaining vis- 
cous solutions for both configurations using the General 
Aerodynamic Simulation Program (GASP), version 2.0. 6 
GASP version 2.0 can solve the full Reynolds-averaged 
Navier-Stokes equations as well as subsets of these equations, 
including the thin-layer Navier-Stokes (TLNS), parabolized 
Navier-Stokes (PNS) and Euler equations. The solutions 
shown in this study were obtained by solving the TLNS equa- 
tions in the nose region in order to capture the region of sub- 
sonic flow behind the bow shock and the PNS equations over 
the remaining configuration. Solutions were obtained for the 
waverider and flat-top configurations over an angle-of-attack 
range from 0° to 10° at Mach 4.0. Off-design solutions at 
Mach 3.5 and Mach 4.5 at 0° angle of attack were obtained 
for the waverider only. 
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Figure 2. Mach 4 Waverider Model 
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Figure 3. Mach 4 Reference Flat-Top Model 
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The numerical method in GASP is based on the 
upwind/relaxation algorithms. Time integration is based on 
the integration of primitive variables and uses a 2-factor 
approximate factorization scheme. Convergence to a steady 
state solution is accomplished by reducing the L2 norm of the 
residual by 4 orders of magnitude. Van Leer’s flux-vector 
splitting algorithm is used with the exception that full flux is 
enforced in the marching direction for PNS solutions. Mesh 
sequencing is used to accelerate convergence for the TLNS 
solutions. A no-slip boundary condition with a fixed wall 
temperature of 324.67 K is imposed on solid boundaries. 

The grid used for the waverider solutions was gen- 
erated from a previous computational study which examined 
this configuration at the design Mach number only. 7 The grid 
was generated using an algebraic transfinite interpolation 
method and used an adaptive grid approach to cluster cells in 
the region of shock waves at the design point. The grid that 
was adapted for flow at Mach 4.0 and an angle of attack of 0° 
is used for all of the solutions in this study. The waverider 
grid consisted of 53 points in the streamwise direction, 71 
points in the circumferential direction, and 51 points in the 
vertical direction. The grid for the flat-top solutions was gen- 
erated using algebraic transfinite interpolation methods with 
elliptic interior point refinement. These techniques were 
implemented by using the GRIDGEN software package. 8 
The flat-top grid contained 51 points in the streamwise direc- 


tion, 71 points in the circumferential direction, and 61 points 
in the vertical direction. Both grids had points clustered near 
the solid boundaries in order to adequately resolve boundary 
layers in viscous solutions. The amount of clustering 
required is determined by examining the inner law variable, 
y + . Previous studies have shown that y + values on the order 
of 1.0 will yield an accurate solution. 9 A diagram of both the 
waverider and flat-top model grids is shown in figure 4. 

Lift and drag coefficients were obtained by integrat- 
ing pressures predicted by CFD solutions over the configura- 
tion surfaces and including an estimate for skin friction using 
the reference temperature method. 10 Additionally, two meth- 
ods are used to account for the blunt base. The first method 
assumes a base pressure of zero while the second assumes a 
base pressure equal to freestream pressure. Lift and drag pre- 
dictions are presented using both methods. 

Results 

The results from the computational study are pre- 
sented in two parts. First, lift and drag coefficients from inte- 
grated surface pressure predictions are presented for the 
design Mach number of 4.0 and selected off-design Mach 
numbers. Second, flow-field solutions as well as lift and drag 
coefficients from integrated surface pressure predictions are 
presented for both the flat-top model and the waverider at 



Figure 4. Grids for Waverider and Flat-Top Configurations 
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Mach 4.0 and various angles of attack. Comparisons of sur- 
face pressure predictions and performance characteristics of 
each configuration are made. Two different methods of 
accounting for the blunt base is also examined for both mod- 
els. Comparisons between computational predictions and 
experimental data are presented for some cases. The experi- 
mental data were documented assuming that the base pres- 
sure is equal to freestream pressure. Therefore, comparisons 
are presented only for this case. 

Lift and drag coefficients for the waverider model at 
0° angle of attack for the design Mach number of 4.0 and for 
off-design Mach numbers are shown in table 1 . The predic- 
tions for lift coefficient, drag coefficient and lift-to-drag ratio 
are presented in figures 5, 6 and, 7. These results show that 
lift decreases as Mach number increases. Drag also decreases 
as Mach number increases for the case of base pressure equal 
to zero. However, for the case of base pressure equal to 
freestream pressure, the drag has a maximum value at Mach 
4.0. Because of this difference, there is an opposite trend in 
the lift-to-drag ratios when comparing the two methods of 
base pressure correction. Assuming zero base pressure, the 



Figure 5. Mach number effect on C L for waverider from 
CFD Predictions at oc=0.0. 



Figure 6. Mach number effect on C D for waverider from 
CFD Predictions at a=0.0. 
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Figure 7. Mach number effect on L/D for waverider from 
CFD predictions at oc=0.0. 


lift-to-drag ratio increases as Mach number increases. How- 
ever, the lift-to-drag ratio decreases with increasing Mach 
number when freestream pressure at the base is assumed. In 
both cases, the lift-to-drag ratios presented agree with refer- 
ence 3, which showed that the performance of waveriders 
does not degrade significantly at off-design conditions. 



Phase = 0 

P — P 

r base r oo 

M 

OO 

C L 

c D 

L/D 

C L2 

C D2 

L/D 

3.5 

0.0823 

0.0218 

3.77 

0.0823 

0.0121 

6.80 

4.0 

0.0747 

0.0189 

3.95 

0.0747 

0.0144 

6.47 

4.5 

0.0677 

0.0165 

4.10 

0.0677 

0.0106 

6.39 


Table 1: Lift and drag coefficients from integrated surface pressure predictions for waverider model for Mach 

4.0 and off-design Mach numbers at 0° angle of attack. 
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A comparison of the waverider and reference con- 
figuration flow fields shows some of the differences between 
the two models. Figure 8 shows solutions of the waverider 
and flat-top configurations at Mach 4.0 and 0° angle of attack. 
Static pressure contours are shown for the centerline and the 
base of each model. The most notable difference is that the 
bottom-surface pressure values are considerably lower for the 
flat-top configuration than for the waverider model. An 
examination of the geometry reveals that the bottom surface 
of the flat-top configuration has a slight expansion, due to the 
method which was used to design this configuration. In con- 
trast, the bottom surface of the waverider model has a com- 
pression surface. The solutions at the base of each model 
show that the shock is slightly detached for each configura- 
tion. This detachment is due to boundary-layer displacement 
and blunt leading-edge effects. The bottom-surface flow field 
is more uniform for the waverider as compared to the refer- 
ence model and thus more attractive for propulsion integra- 
tion. 

An examination of solutions at higher angles of 
attack at Mach 4.0 show the same type of behavior. For 
example, static pressure contours at the centerline and base of 


the reference and waverider models at Mach 4.0 and 4° angle 
of attack are shown in figure 9. The solutions at 0° and 4° 
angle of attack are representative of all the angles of attack 
studied in that for each case the bottom-surface pressure val- 
ues are lower for the flat-top model than for the waverider. 
The shock detachment distance at the leading edge also 
increases as angle of attack increases. 

The lower pressure values on the bottom surface of 
the flat-top configuration result in lower lift and drag values. 
The lift and drag coefficients from integrated surface pressure 
predictions for the waverider and flat-top models at Mach 4.0 
are shown in tables 2 and 3, respectively. The predicted val- 
ues are again presented with two methods of base-pressure 
correction. At all angles of attack studied, the predicted val- 
ues of both lift and drag coefficients are lower for the flat-top 
model than for the waverider. The performance characteris- 
tics of each configuration are examined by presenting the lift 
and drag predictions as a function of angle of attack and the 
lift-to-drag ratios as a function of lift coefficient. 

The predicted and experimental lift and drag charac- 
teristics of each configuration at various angles of attack are 
shown in figures 10 through 14 for the design Mach number 


P/P Contours at Center Line 
Waverider 


P/P m Contours at Base 
Waverider 



Figure 8. TLNS/PNS Solutions at M =4.0, a=0.0° 
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171* Contours at Center Line 
Waverider 


I’/P Contours at Center Line 
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Phase = 0 

P — P 

r base r oo 

a 

C L 

c D 

L/D 

C L2 

C D2 

L/D 

0.0 

0.0747 

0.0189 

3.95 

0.0747 

0.0115 

6.45 

1.0 

0.0933 

0.0218 

4.28 

0.0934 

0.0144 

6.47 

2.0 

0.1136 

0.0264 

4.30 

0.1139 

0.0190 

5.99 

3.0 

0.1344 

0.0307 

4.38 

0.1283 

0.0233 

5.56 

4.0 

0.1553 

0.0364 

4.27 

0.1558 

0.0290 

5.37 

6.0 

0.2196 

0.0554 

3.96 

0.2204 

0.0480 

4.59 

10.0 

0.3155 

0.0982 

3.21 

0.3168 

0.0909 

3.39 


Table 2: Lift and drag coefficients from integrated surface pressure predictions for waverider model at Mach 

4.0 and selected angles of attack. 


of 4.0. The variation of lift coefficient with angle of attack is 
presented in figure 10. This figure shows that the lift coeffi- 
cient values for the waverider model are higher than those for 
the flat-top model at each angle of attack. These results are 
expected based on the predicted pressure values in the com- 


putational solutions. Good agreement is generally obtained 
between the CFD solutions and the experimental values. The 
drag coefficient data are presented in figure 1 1 assuming both 
zero base pressure and freestream base pressure. In both 
cases, higher drag values are shown for the waverider model 
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Phase = 0 

P — P 

r base r oo 

a 

C L 

c D 

L/D 

C L2 

C D2 

L/D 

0.0 

0.0476 

0.0155 

3.07 

0.0476 

0.0081 

5.87 

1.0 

0.0616 

0.0165 

3.73 

0.0617 

0.0092 

6.74 

2.0 

0.0804 

0.0192 

4.19 

0.0807 

0.0118 

6.82 

3.0 

0.0992 

0.0226 

4.39 

0.0996 

0.0152 

6.53 

4.0 

0.1178 

0.0268 

4.40 

0.1183 

0.0194 

6.11 

6.0 

0.1684 

0.0399 

4.22 

0.1692 

0.0325 

5.20 

10.0 

0.2473 

0.0716 

3.45 

0.2487 

0.0644 

3.86 


Table 3: Lift and drag coefficients from integrated surface pressure predictions for Flat-Top Reference model 

at Mach 4.0 and selected angles of attack. 



Angle of Attack (Degrees) 

Figure 10. CFD and experimental C L vs. a at M =4.0 
for waverider and flat-top model. 



Figure 11. CFD Predictions for C D vs. a at M =4.0 


at each angle of attack. The drag values for the waverider 
also increase more rapidly as the angle of attack increases. A 
comparison of drag coefficient values between CFD and 
experimental data for the case of base pressure equal to 
freestream pressure is shown in figure 12. Again, good agree- 
ment is generally obtained between CFD results and experi- 
mental data. 



0.0Q 

0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 

Angle of Attack (Degrees) 

Figure 12. CFD and experimental C D vs. a 
at M =4.0 with P. =P . 

°° base °° 

The lift-to-drag ratio as a function of angle of attack 
at Mach 4.0 assuming zero base pressure is presented in fig- 
ure 13. For this case, the maximum lift-to-drag ratio for the 
waverider occurs at an angle of attack of 3° and for the flat- 
top configuration, at 4.0 degrees. The maximum lift-to-drag 
ratio of the flat-top model is less than 1 percent higher than 
that for the waverider model. The lift-to-drag ratio for both 
configurations are presented in figure 14 assuming a base 





Angle of Attack (Degrees) 

Figure 13. L/D vs. a for waverider and flat-top model 
with P base =0.0 from CFD predictions. 


pressure equal to freestream pressure. Using this method, the 
predicted maximum lift-to-drag ratio occurs at 1° angle of 
attack for the waverider and at 2° angle of attack for the flat- 
top model. The predicted maximum lift-to-drag ratio for the 
flat-top model is approximately 5.4 percent higher than that 
of the waverider model. The experimental values for lift-to- 
drag ratio are also shown in figure 14. In general, the lift-to- 
drag values from CFD solutions are higher than the experi- 



Angle of Attack (Degrees) 

Figure 14. Experimental and CFD predictions for 
L/D vs. a at M =4.0 with P. =P . 

oo base °o 


ment, especially for the flat-top model and the same trends are 
observed in the experiment and CFD results. 

The lift-to-drag ratio as a function of lift coefficient 
is presented in figure 15 assuming a base pressure equal to 
zero and in figure 16 assuming a base pressure equal to 
freestream pressure. The experimental data are included in 
figure 16 for comparison purposes. For small values of lift, 
the flat-top model has a higher lift-to-drag ratio than the 
waverider. For lift coefficient above 0.17, the waverider 
shows better performance. These results also show that the 
flat-top model has a higher maximum lift-to-drag ratio, but it 



with P base =0.0 from CFD predictions. 



C L vs. L/D with P base =P„. 
occurs at a lower lift coefficient. 

Interpretations 

Several conclusions can be made by the comparison 
between computational solutions of both configurations and 
between computational solutions and experimental data. The 
waverider configuration has higher lift than the reference flat- 
top configuration for all conditions investigated. However, 
the flat-top model has a higher lift-to-drag ratio at some 
points and a higher maximum lift-to-drag ratio than the 
waverider model due to relatively lower drag. In other words, 
at most values of lift coefficient, the drag coefficient is lower 
for the flat-top model than for the waverider, due to the lower 
bottom-surface pressures. This results in a higher lift-to-drag 
ratio for the flat-top model at most conditions. The improved 
performance of the flat-top model over that of the waverider 
is not due to the flat top surface, but rather to the shape of the 
bottom surface. For the configurations investigated herein, 
the bottom surface of the flat-top model has a slight expan- 
sion and therefore, lower surface pressure. 
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The results also illustrate that the type of base pres- 
sure assumed has a significant effect on the relative perfor- 
mance of the two configurations. The differences in lift-to- 
drag ratio are much greater when freestream pressure is 
assumed at the base than when zero base pressure is assumed. 
In comparing two configurations to assess their suitability for 
integration into a practical hypersonic vehicle, the assump- 
tion of zero base pressure may be more appropriate since the 
blunt base will most likely be eliminated in any representative 
vehicle design by the addition of control surfaces. However, 
it would be appropriate to use the assumption of freestream 
pressure acting at the base in the design process in order to 
avoid widely varying solutions due to changing base areas. 

The flat-top model exhibits the same shock attach- 
ment properties as the waverider model because the planform 
shapes of the two configurations are identical. Based on the 
improved performance of the flat-top model, it can be con- 
cluded that the basic physical effect that gives the waverider 
its advantage in aerodynamic performance over conventional 
hypersonic vehicles is the shock attachment caused by the 
shape of the leading edge. Designing the bottom surface as a 
streamsurface improves flow-field uniformity and makes the 
waverider a more attractive candidate for scramjet integra- 
tion, but does not enhance the aerodynamic performance. 
The bottom surface of a waverider can be altered somewhat 
without a large decrease in lift-to-drag ratio. This suggests 
that altering the lower surface of a waverider configuration to 
integrate engine components should not degrade the perfor- 
mance of the vehicle significantly as long as the planform 
shape is maintained. 

Conclusions 

A computational study was conducted to better 
understand the results obtained from an experimental investi- 
gation of a Mach 4 waverider and a reference flat-top config- 
uration. The results from this experiment contradicted 
waverider design theory by showing that the flat-top configu- 
ration had a higher maximum lift-to-drag ratio than the 
waverider. The waverider model was designed using an opti- 
mization routine which includes an estimate for skin friction 
in the optimization process. The reference model was 
designed by adjusting the top surface of the waverider to 
obtain a flat surface, while maintaining the same planform 
shape and cross-sectional area distribution. 

Viscous solutions were obtained for both configura- 
tions and lift and drag predictions were obtained from these 
solutions by integrating surface pressure values and including 
an estimate for skin friction. Solutions of the waverider at the 
design Mach number and selected off-design Mach numbers 
showed that the off-design performance of the waverider does 


not degrade significantly at off-design conditions. The bot- 
tom-surface pressures of both the waverider and flat-top 
model at Mach 4.0 and selected angles of attack showed that 
the bottom surface of the flat-top model provides a slight 
expansion, in contrast to the waverider bottom surface, which 
acts as a compression surface. The result is that the lift and 
drag predictions for the flat-top model are much lower than 
those for the waverider model. These characteristics result in 
higher lift-to-drag ratios for the flat-top model for lift coeffi- 
cients lower than 0.17 and a higher maximum lift-to-drag 
ratio. However, the differences in performance between the 
two configurations are primarily due to the shape of the bot- 
tom surface independent of the top surface. The results also 
indicate that the performance of each configuration is signif- 
icantly affected by the type of base pressure correction used. 

The CFD results also show that the reference config- 
uration exhibits the same shock attachment properties as the 
waverider, due to the planform shapes being identical. The 
results from this study suggest that the leading-edge shock 
attachment, which is a result of the planform shape, is the 
main effect that gives waveriders their high lift-to-drag ratios. 
The performance of the configuration should not degrade sig- 
nificantly if the bottom surface is altered from the optimized 
waverider design, as long as the planform shape is main- 
tained. 
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